Evidence for electron-electron interaction in topological insulator thin films 
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We consider in our work single crystal thin films of Bi2Se3, grown by molecular beam epitaxy, both with 
and without Pb doping. Angle-resolved photoemission data demonstrate topological surface states with a Fermi 
level lying inside the bulk band gap in the Pb doped films. Transport data show weak localization behavior, 
as expected for a thin film in the two-dimensional limit (when the thickness is smaller than the inelastic mean 
free path), but a detailed analysis within the standard theoretical framework of diffusive transport shows that 
the temperature and magnetic field dependences of resistance cannot be reconciled in a theory that neglects 
inter-electron interactions. We demonstrate that an excellent account of quantum corrections to conductivity is 
achieved when both disorder and interaction are taken into account. These results clearly demonstrate that it 
is crucial to include electron electron interaction for a comprehensive understanding of diffusive transport in 
topological insulators. While both the ordinary bulk and the topological surface states presumably participate in 
transport, our analysis does not allow a clear separation of the two contributions. 



I. INTRODUCTION 

Topological structures in condensed matter often result 
from interactions; examples include vortices in superconduc- 
tors, superfluids, two-dimensional XY systems, and the frac- 
tional quantum Hall effect. Topology also governs the be- 
havior of certain non-interacting systems: the integer quan- 
tum Hall effect is the quintessential example, while the re- 
cently proposed topological band insulators provide another. 
The latter are predicted to contain conducting surface states 
with an odd number of Dirac cones that are topologically pro- 
tected against time reversal invariant perturbation. It is im- 
portant to demonstrate that the topological features of these 
new systems are robust to interparticle interaction. The exis- 
tence of surface states has been investigated and confirmed for 
several candidate topological insulators using angle-resolved 
photoemission spectroscopy (ARPES).^'^ Detecting these sur- 
face states by transport, which might at first appear the most 
natural probe, has proved more challenging, largely because 
of non-negligible bulk conduction. Here, we report a detailed 
study of the temperature and magnetic field dependence of 
the conductivity in topological insulator thin films grown by 
molecular beam epitaxy. A theoretical analysis of the ob- 
served quantum corrections to the conductivity unambigu- 
ously demonstrates that electron-electron interactions play a 
crucial role in determining the diffusive transport behavior of 
topological insulators. 

Bi2Se3 is a candidate for a three dimensional topological 
insulator (TI), provided its Fermi level can be shifted to lie 
inside the bulk energy gap, as can be accomplished by appro- 
priate amount of doping. Due to its relatively large band gap 
(0.3 eV) and its simple surface state structure (a single Dirac 
cone), it has become the canonical reference material for the 
3D TIs.^^ Several transport experiments have suggested the 
presence of conducting surface states in candidate 3D TIs. A 
large magnetofingerprint signal^ in macroscopic crystals of 




FIG. 1. A Scanning Tunneling Microscope (STM) Image of a Typ- 
ically Grown Bi2Se3 Film with a Thickness of 45 Quintuple Layers 
(QLs). 



nonmetallic Bi2Se3, Aharonov-Bohm oscillations in Bi2Se3 
nanoribbons"^ with dimensions of ^ 100 nm, and Shubnikov- 
de Haas (SdH) oscillations" in bulk TI (Bi2Te3) have been 
reported and associated with the surface states. However, an- 
other SdH measurement on high quality Bi2Se3 single crys- 
tals did not find significant contributions from surface states 
and suggested that these are subject to strong scattering.'^ It 
is noted that the analysis of all the above experiments neglects 
interactions. 

Recent progress in thin film growth of TI materials by 
molecular beam epitaxy (MBE)'-' '"' has enabled the fabrica- 
tion of thin films of crystalline Bi2Se3. We study the tem- 
perature and magnetic field dependence of two thin films of 
Bi2Se3, one undoped and the other doped. The latter is doped 
with an appropriate concentration of Pb, which should ideally 
render the bulk insulating. ARPES measurements demon- 
strate that the latter thin film possesses Dirac-cone surface 
states with the Fermi level lying in the bulk band gap. Al- 
though it is likely that there is bulk conduction in the doped 
thin film as well, the expectation is that the contribution from 
surface states will make a significant contribution to the elec- 
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FIG. 2. ARPES spectra of the 45 nm thick crystalUne films, (a) The 
undoped Bi2Se3 film and (b) doped Bi2-iPba;Se3 film. The energies 
are measured relative to the Fermi level. 



trical transport. The conductance shows a logarithmic tem- 
perature dependence at low temperatures, a hallmark of weak 
localization (WL) in two dimensions (2D), originating from 
the interplay between disorder, spin-orbit coupling and/or 
electron-electron interactions (EEI)i^. We find that the tem- 
perature and magnetic field dependences of our data are qual- 
itatively inconsistent with conventional 2D transport theory of 
non-interacting electrons, but are successfully explained when 
we include, following the theory of Lee and Ramakrishnan'^, 
quantum corrections to the conductivity originating from EEI 
and also account for the Zeeman splitting, which is large in 
these narrow band gap materialsii. 



II. EXPERIMENTAL RESULTS 

We use MBE to grow 45 QL Bi2Se3 films on bare insulating 
6H-SiC (0001) substrates whose resistivity is ^ 1 x 10^ Q-cm, 
allowing us to neglect the conductance contribution from the 
substrate. A scanning tunneling microscope (STM) image of 
a typical film is shown in Figure [T] Atomically flat terraces 
with widths over 100 nm can be observed, with 1 quintuple 
layer thick terraces. Figure|2|a) shows the ARPES band map 
of a 45 QL Bi2Se3 film. A single Dirac cone is observed at 
the r point, with the Dirac point located at 0. 1 85 eV below the 
Fermi level. Hence, both the surface states and the bulk con- 
duction band can contribute to electronic transport. The bulk 
chemical potential can be tuned inside the bulk energy gap 
by using appropriate dopants.^ In this experiment, we accom- 
plish this by doping the sample with .37% Pb.'^ As shown 
in Fig.|2b), after doping, the Fermi level of the Bi2-KPb2;Se3 
film is inside the energy gap. In this situation, the bulk con- 
ductivity is suppressed and the surface conductance becomes 
more evident. 

The Bi2Se3 films were grown under Se-rich conditions on 
6H-SiC (0001) substrates at 220°C in an ultrahigh-vacuum 
(UHV) system (Omicron), equipped with MBE, STM and 
ARPESii. The base pressure of the system is 1.5 x 10^^" 
Torr. High purity Bi (99.9999%) and Se (99.999%) were ther- 
mally evaporated from standard Knudsen cells. The tempera- 
tures of the Bi source and Se source are 550 and 170°C. The 
Se4(Se2)/Bi flux ratio was between 10 and 15, which leads 
to a growth rate of ^ 0.3 QL min^^. In ARPES measure- 
ment. He-la (21.21 eV) photons produced by a Gammadata 
VUV 5000 discharging lamp and a Scienta SES2002 analyzer 
are used to excite and collect photoelectrons, respectively. All 
STM and ARPES data are taken at room temperature. 

Subsequent to the ARPES analysis, the films were covered 
by 30 nm thick amorphous Se, which is quite insulating, as 
protective layers. Then, the samples were taken out from the 
UHV system for transport measurements in a variable tem- 
perature magnetocryostat with high magnetic field (H < 80 
kOe) and low temperature (0.1 K < T < 300 K) capability. 
Standard four probe electrical transport measurements were 
carried out using an ac resistance bridge with an excitation 
current of 500 nA. All resistance values were obtained by av- 
eraging over 50 measurements. 

We determine the carrier density in the two samples using 
the Hall effect, as shown in Fig. [3] The carriers are electrons 
in both samples, and from the approximately linear behavior 
of Rxy at low magnetic fields, we estimate a carrier density of 
2.27 X 10^^ cm-3 (1.02 x 10" cm^^) for the BizSeg film and 
1.1 X 10^^ cm-3 (4.95 x lO^^ cm^^) for the Bi2-2:Pba;Se3 
film at 2 K. We estimate kpie ^ 22 for Bi2Se3 and 5.7 
for Bi2_a:Pba;Se3, well within the diffusive transport regime. 
(The quantity £e is the elastic scattering length.) R^y has a 
slightly non-linear behavior at high magnetic fields, which ap- 
pears to be quite generic at low temperatures, and is likely an 
indication of more than one transport channel. '^'^^ 

The temperature dependence of the sample resistivity in the 
range 0.1 K < T < 300 K reveals a qualitative difference be- 
tween the electrical properties of the two thin films. See Fig.|4] 
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FIG. 3. Hall Resistance of the 45 QL Bi2Se3 and Bi2-a:PbxSe3 
Films at 2 K. Below 20 kOe, the Hall resistance shows almost linear 
behavior (the red solid lines in the figure) in both samples. However, 
with the increasing magnetic field, we can clearly see the deviation 
from the linear behavior, which is larger in the Bi2-a:Pba;Se3 film. 
The nonlinear Hall behavior is reminiscent of the two channels of 
carriers in the films. One channel may be from the surface state. 



for the resistance vs temperature (R-T) plots for the two films 
from room temperature to 0.5 K. Above 50 K, the resistance 
of the Bi2Se3 film increases linearly with temperature, char- 
acteristic of the metallic behavior expected for a degenerately 
doped semiconductor Its resistivity varies from 1 .43 mfi cm 
at 297 K to 0.72 mQ cm at 500 mK. For the Bi2-:EPbj,Se3 film, 
the resistance is higher, with only a very weak, non-monotonic 
T dependence. Its resistivity is 3.69 mil cm at both 297 K and 
500 mK. 

In an earlier studyi^ of a bulk Bi2Se3 system, it was found 
that the bulk remains conductive even for densities much 
smaller than those quoted above. The origin of the qualita- 
tive difference between our two thin films is not understood, 
but presumably has to do with the low-dimensional nature of 
our system, and also the fact that it has a much lower mobility 
than the system in Ref j^^. 

We have also measured the magnetoresistance for fields 
up to 80 kOe for three orientations of the magnetic field, as 
shown in Fig. |5] In the perpendicular field {H±), the frac- 
tional MR change, AR/Rq = (i?(i?) - R{0))/Ro, between 
80 and kOe is about 35% for the BiaSea film and 25% for 
the Bi2-2;Pbj;Se3 film. When H is parallel to the film, the frac- 
tional MR change is much smaller, between 2.7% and 5.9%. 
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FIG. 4. Resistance vs Temperature Plots of the 45 QL Bi2Se3 and 
Bi2-:cPb^Se3 Films from 300 K to 0.5 K. 



Linear MR was observed above 10 kOe (30 kOe) when a 
perpendicular field was applied to the Bia-^PbrSea (BiaSea) 
film as shown in Fig. ISjb). Such a linear MR has been at- 
tributed to the quantum linear MR of the surface states^. We 
note that even at room temperature, the MR above 45 kOe 
has a better fit with a linear rather than parabolic dependence 
on the magnetic field for both samples (see Fig. |6ll. Upon 
closer inspection, the MR of the BiaSea and Bi2-:EPb2;Se3 
films shows sharp minima at zero field at low temperatures. 
These are better revealed in the Figure [T] 

In the following, we focus on the behavior at low temper- 
atures and low magnetic fields, because that is the most rele- 
vant parameter regime for our discussion within the context of 
quantum theories of localization. Figures [HI a) and (b) display 
the T dependence of the conductance for zero magnetic field 
as well as for a 20 kOe in-plane field. The logarithmic be- 
havior is suggestive of a weakly disordered 2D system where 
diffusive transport is determined by either weak localization 
or EEL Figures |8jc) and[8jd) show the magnetoconductance 
as a function of a perpendicular field H±, while Figs. |8je) 
and [HI f) show the magnetoconductance for the two directions 
of an in-plane field (i7|| and H'^) as defined in the caption of 
Fig. [5] The observed magnetoconductance is negative, which 
is often referred to as weak anf/localization. In all cases, the 
magnetoconductance peak disappears at temperatures above 
20 K, suggesting quantum correction as its origin. 
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FIG. 5. Fractional change in magnetoresistance, AR{H) = 
{R{H) - J?(0))/i?(0), of the 45 nm thick (a) undoped Bi2Se3 and 
(b) doped Bi2-i:Pbi:Se3 films. We used three different orientations 
of the magnetic field: denotes magnetic field perpendicular to 
the surface of the thin film, while Hw and iJy denote an in plane 
magnetic field perpendicular to and parallel to the excitation current, 
respectively. For the latter two, the resistance ratio has been multi- 
plied by the factor shown on the plot. All data are at 500 mK. 



III. ANALYSIS 

A. Weak localization for noninteracting electrons 

We first attempt to analyze these measurements using the 
standard results of weak localization (WL) theory for nonin- 
teracting electrons, and find that they are inconsistent with ei- 
ther 3D or 2D weak localization. The logarithmic temperature 
dependence of conductance rules out 3D behavior, because in 
3D we expect Actwl oc T~p/^, where p is defined by the T 
dependence of the phase coherence time cx T"^'. The ex- 
ponent p is dependent on the source of the inelastic scattering 
which causes the phase decoherence. For example, -p — Z for 
phonon scattering, and -p — Iji (3D) or p = 1 (2D) for the 
electron-electron interaction. In 2D, on the other hand, the 
correction to the conductivity at zero field is given by2^ 



AawL(T,i/ = 0) 



(1) 



where a is a constant depending on the relative strengths of 
the spin-orbit and spin-flip (magnetic) scatterings. In the limit 




-80-60-40-20 20 40 60 80 
H (kOe) 

(b) 770r 




-80-60-40-20 O 20 40 60 80 
H (kOe) 

FIG. 6. Resistance vs Magnetic Field Plots at 300 K for the 45 QL 
Bi2Se3 and Bi2-iPba;Se3 Films in the Perpendicular Field. 



of weak spin orbit and magnetic scattering, one obtains a = 1; 
in the limit of strong spin orbit scattering and weak magnetic 
scattering, one finds a — —l/'i,; and when the magnetic scat- 
tering is strong, we have a — Q. The Hikami-Larkin-Nagaoka 
equation for magnetoresistance in perpendicular field is^^. 



A(7wL(^f±) - A(7wl(0) 

X -0 



he 



In 
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where ip is the digamma function and £^ is the inelastic scat- 
tering length. The positive slope of conductance vs. InT 
implies that our samples are in the limit of weak spin-orbit 
scattering (a = 1). However, positive a implies positive mag- 
netoconductance, which is qualitatively inconsistent with the 
observed anti-localization behavior. 

Figure|7]shows explicit comparison of the experimental re- 
sults to the weak localization theory. The lines in Figure |7] 
(c-d) show the weak localization theory (Eq. (2) in the main 
text) for a = 1, as implied by the positive slopes of the data 
in Figure|2](a-b). We have taken a typical value for the coher- 
ence length, £^ — 500 nm. The lines in Figure |7](e-f) are fits 
to the weak localization theory of Maekawa and Fukuyamal^ 
in parallel field. 
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FIG. 7. Weak localization theory, (a-b) Temperature dependence of 
the conductance. Lines are fits to the weak localization theory in 2D 
(solid) and 3D (dashed and dash-dotted, corresponding to phonon 
and electron-electron scattering as the source of dephasing). The 
temperature dependence of both films suggests that our films are in 
the regime of weak spin-orbit scattering. However, in this regime, 
the theory predicts a positive magnetoconductance in perpendicu- 
lar field, (c-d) Change in conductance, AG = AG(//) - AG(0). 
The theory predicts a localization effect, however we observe anti- 
localization, (e-f) Magnetoconductance in parallel fields {H\\ de- 
notes a field parallel to the film and perpendicular to the excitation 
current, while //y denotes a field parallel to both the film and the cur- 
rent). Lines are best fits to the theory of Maekawa and Fukuyama.^^ 
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g is the Zeeman g-factor, tq refers to the scattering time from 
disorder, and Tso,i refers to the spin-orbit scattering in the di- 
rections perpendicular and parallel to the film. The source 
of the magnetic field dependence is the Zeeman splitting only, 
and this effect vanishes for 5^-0. Roughly speaking, the fac- 
tor {y + z)/{x ~ y) is a measure of the comparative strengths 
of the spin orbit and elastic scattering. For example, if elastic 
scattering is strong (large x) compared to spin-orbit scattering, 
then this factor will be small. In the low-field limit, this equa- 
tion gives Act cx B^, which is different from the observed 



sharp cusp. Nevertheless, as shown in Fig.|7te) and (f), a very 
narrow parabola can be obtained for [y + z)/{x~ y) « 0.01; 
i.e. weak spin orbit scattering {x ~ y ^ y ^ z). Although 
consistent with the temperature dependence, it is inconsistent 
with the expectation that the spin-orbit interaction in Bi2Se3 
is strong. 



B. Electron-electron interaction 

Having shown that spin-orbit related WL alone cannot ex- 
plain our data, we now examine the role of the dynamically 
screened EEI which also produces corrections to the conduc- 
tivity. Again the logarithmic T dependence rules out the 3D 
theory where we expec t^^'^^ AaEEi ^/T. On the other 
hand, the correction in 2D is 



AaEEi(T) 



In ( — 

To 



(8) 



where Fo- is a function of the average of the static screened 
Coulomb interaction over the Fermi surface, F = J dCl v{q = 
2kFsm{e/2))/ JdClv{0). In2D, f'a = -4+(8/F)(l + 
F/2)lii(l + F/2). A non-zero magnetic field activates the 
large Zeeman splitting in Bi2Se3^, introducing a further cor- 
rection to the conductivityji^ 



AaEEi(T, if) - AaEEi(T, 0) 



where 



92{T,H) 



dn In 



F,g2{T,H), (9) 



dw- 



- 1 

(10) 

In Bi2Se3, the Zeeman g-factor is 23 when H is in-plane, and 
32 when H is out of plane'^. When the Zeeman energy is 
much larger than fcsT, 32 (i?) ~ Mg^isH / {l.'ikBT)). At 
20 kOe (red data shown in Figs. jOfa) and (b)), this limit is 
satisfied for T <C 3 K. The conductance is related to con- 
ductivity by a geometrical factor (e.g. length/width for a 2D 
surface) ^1 for our samples. 

Figure [8] shows a comparison of the data with the electron- 
electron interaction (EEI) theory. There are two fitting param- 
eters, which are both obtained from the slopes of the lines in 
Figure [8] (a-b). These are the screening parameter, F^, and a 
geometrical factor relating the conductance with the conduc- 
tivity, G = fa. In a wire with uniform current, / is the ratio 
of cross-sectional area to length. We do not make this assump- 
tion, and instead use this factor as a parameter of the model. 
Using no further parameters, the EEI expressions correctly re- 
produce the signs of the T and H dependences, which already 
indicates the importance of EEI interactions. Although it does 
not capture the sharp peaks at zero field, the theory nicely re- 
produces all of the other qualitative features. 
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FIG. 8. Electron-electron interactions theory, (a-b) Temperature de- 
pendence of the conductance. Lines are fits to the 2D (solid) and 3D 
(dashed) theory of Lee and Ramakrishnan._ The slope of the solid 
lines determines the screening parameter (Fa = 0.47 and 0.75 for 
undoped and doped samples). Change in magnetoconductance in (c- 
d) a perpendicular magnetic field and (e-f) parallel fields. Solid lines 
in (c-f) are fits to the electron-electron interaction theory using no 
additional fitting parameters besides those found in (a-b). 



C. Combined weak localization and electron-electron 
interaction 

We finally proceed to fit the experiment to the theory that 
contains corrections due to both weak localization and inter- 
action. The temperature and dependences in Figure |9ja)- 
(d) are fitted to the combined model, using the approximation 
Aa{T, H^) = A(TEEi(r, Hx_) + A(7wL(r, i/^), with F,, a 
and ^(f, taken as fitting parameters. (We take p — 1, as appro- 
priate for 2D systems at sufficiently low temperatures where 
EEI is the dominant mechanism for dephasing.) The data were 
fitted in the range —2 kOe< H <2 kOe, which corresponds to 
the "low-field" region of the Bi2Se3 plot. We find a = —0.31 
(-0.35), = 1.1 (0.64) ^m, and F„ = 0.15 (0.35) for 
Bi2Se3 (Bi2-2:Pb2;Se3). The values for a are slightly different 
from the value of — 1 /2 expected from the strong spin-orbit in- 
teraction. The inelastic mean free paths are consistent with a 
previous measurement, = 0.5/im at 2 Ki^'. 

The parameter F^, which incorporates the effect of in- 
teraction through F, can be estimated from first principles. 
While our sample thickness is much smaller than the inelas- 
tic scattering length, making it two dimensional insofar as the 
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FIG. 9. Transport properties in low temperature and low magnetic 
field. The left columns are for the undoped Bi2Se3 film, while the 
doped Bi2-a;Pba;Se3 data is on the right, (a-b) Conductance vs T for 
zero magnetic field (black circles) and for = 20 kOe (red trian- 
gles), with fits to 2D theories (solid lines), and 3D electron-electron 
interaction (EEI) theory (dash-dotted lines), (c-d) Magnetoconduc- 
tance AG = G{H) — G(0) for fields perpendicular to the film. Solid 
lines are the result of a combined WL and EEI theory, as explained in 
the text, (e-f) Magnetoconductances for in plane fields H'^ and i/y , 
parallel and perpendicular to the current direction. The solid lines are 
fits to the combined WL and EEI theory. The magnetoconductance 
data are all taken at 500 mK. 



weak localization correction to conductivity is concerned^^, 
the thickness is much larger than the Thomas Fermi screening 
length (fcxF < 1 nm), suggesting that F is to be estimated 
in a three-dimensional calculation. We will assume that to be 
the case. For the carrier densities reported above, a straight- 
forward calculation using the Thomas-Fermi approximation 
for the interaction gives F = 0.68 (F„ = 0.62) and 0.72 
{F„ = 0.65) for the BiaSeg and Bi2-:rPbj;Se3 films. For 
the calculation of F we have assumed a three dimensional 
parabolic dispersion relation /2m* where to* — 0.15me 
(TOe is the charge of an electron)^"*. We have used the same 
dispersion for both samples, assuming that the doping has a 
negligible effect on the effective mass. (We note that if we 
had taken 2D screening with parabolic dispersion, we would 
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get for F 0.60 and 0.66, respectively. A 2D calculation as- 
suming a Dirac cone with Fermi velocity vp = 5 x 10^ m/s 
yields a calculated F — 0.77 independent of the Fermi wave 
vector.) The fitted values of F{F) are smaller, suggesting that 
we are overestimating the strength of the interaction. This 
discrepancy is not understood, but is not surprising in view 
of the crudeness of our model; specifically, we neglect: the 
strong anisotropy of the effective mass (which can change by 
a factor of 8 between parallel and perpendicular directions 
doping dependence of the mass (which can change its value 
by as much as a factor of twc^^); and corrections to the mass 
due to quantum confinement in one direction. 

We have also fit the data in parallel field, shown in Fig.|9je- 
f). Since the theory does not allow for a difference between 
the two in-plane field directions, we have performed a simul- 
taneous fit to both data sets. While the two data sets agree 
well at low fields, the disagreement at higher fields is not un- 
derstood. In the limit of large spin-orbit scattering and low 
disorder, the predicted behavior cannot explain the sharp peak 
observed in the experiments. However, when the spin-orbit 
and disorder scattering times are comparable (roughly speak- 
ing, this occurs in the parameter range {x + z)/{x — y) ^ 1), 
the weak localization model (Eqn. (O) gives a good account 
of the observed behavior. 

We note that there have been other experimental studies of 
thin films of TI materials^^^**. These do not identify the es- 
sential role of electron electron interactions. After the comple- 
tion of this work we became aware of another article^S which 
studies transport in much thinner Bi2Se3 films (1-6 QLs) and 



interpret their results as indicating EEI as the origin of local- 
ization; their samples are not doped, however, and they also 
did not fit their data to EEI corrections discussed in the present 
article. 



IV. CONCLUSION 

In conclusion, we have studied two thin film samples of the 
Bi2Se3, one of which has surface states with Fermi energy 
located in the conduction band, while the other has a Fermi 
energy within the bulk band gap because of Pb doping. We 
find a 2D WL theory combined with EEI fits the transport 
data of both samples very well. Note that the 2D WL model 
is relevant because the inelastic length in the these samples is 
much larger than the sample thickness,'^ and does not in itself 
indicate transport in a surface state. Although ARPES results 
indicate that the surface states should play a role in the con- 
ductance of the doped sample, it is not possible at this stage 
to single out the contribution from the surface states. What 
we can conclude with confidence is that the electron-electron 
interaction, which has been neglected in most previous analy- 
ses, plays an essential role in quantum corrections to transport 
in thin films of TI Bi2Se3 material. 
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